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PARAMETRIZATION OF POLARIZED PARTON DISTRIBUTIONS 



S. KUMANO 
Department of Physics, Saga University 
Saga 840-8502, Japan 



We report the polarized parton distributions proposed by the Asymmetry Analysis Collab- 
oration (AAC). Using parametrized distributions at Q 2 =l GeV 2 and measured A\ data, we 
determine optimum polarized distributions in the leading order (LO) and next-to-leading order 
(NLO). We find that the polarized antiquark distribution is not well determined particularly 
at small x. It could lead to a rather small quark spin content in comparison with usually- 
quoted values of 10~30%. In our analysis, it varies from 5% to 28% depending on the small- a; 
extrapolation. It is necessary to have small x (< 10 -3 ) data for precise determination. In ad- 
dition, the large- x region should be also studied for Aq, which cannot be determined solely by 
gi. We propose three sets of distributions as the longitudinally polarized parton distribution 
functions. 



1 Introduction 

From the measurements of the proton's polarized structure function g\, we learned that the 
proton spin cannot be understood in a simple quark model. It should be interpreted mainly by 
a gluon contribution and/or effects of angular momenta; however, the situation is not satisfac- 
tory for specifying the major carrier of the proton spin. Our study is intended to understand 
the current status of polarized parton distributions by analyzing inclusive spin asymmetry A± 
data prior to RHIC-Spin measurements. Semi-inclusive data became available, but they are not 
accurate enough to provide any significant constraint. Our group name AAC stands for Asym- 
metry Analysis Collaboration, which consists of theorists and experimentalists0 In this paper, 
we report our analysis of the asymmetry A\ for getting optimum polarized parton distributions. 

In Sec.^, a parametrization form of the polarized parton distributions is discussed. The 
results are shown in Sec.|3[ According to our analysis, the polarized antiquark distribution 
cannot be well determined, and this issue is discussed in Sec.||. The results in Sec.|| and || are 
quoted from Ref.l. Finally, conclusions are given in Sec.|5[ 



2 Parametrization 



First, an x-dependent functional form of parametrized distributions is explained. Because the 
positivity condition is imposed not only in the leading order (LO) but also in the next-to-leading 
order (NLO), it is convenient to have an unpolarized distribution multiplied by an x-dependent 
function: 

Afi(x, Ql) = Ai x a > (1 + lt x x >) fi(x, Ql), (1) 

at the initial Q 2 point (Qq=1 GeV 2 ). Since there is no accurate data to find the difference 
between Au, Ad, and As, these distributions are assumed to be the same. Of course, they are 
not expected to be equal by considering the unpolarized situation! With this assumption and 
the flavor number Nf = 3, we have four distributions (Au v , Ad v , Aq, Ag) to be determined by a 
X 2 analysis. Because there are four parameters for each distribution, we have sixteen parameters 
in total. However, the first moments of Au v and Ad v are fixed (r] Uv = 0.926 and = —0.341) 
by using semi-leptonic decay constants. It means that the number of actual free parameters is 
fourteen. 

The polarized distributions are evolved! to experimental Q 2 points of the spin asymmetry 
At, which is expressed by the structure functions as 

, . q2) „ 2x 9l {x,Q 2 )[l + R{x,Q 2 )] 

The longitudinal-transverse structure function ratio R is taken from the analysis of SLAC-1990, 
F2 is calculated by the GRV98 distributions, and g\ is calculated by our parametrized distri- 
butions. Then, x 2 1S evaluated in comparison with experimental data: x 2 = XX^i ata ( x > Q 2 ) — 
A1 alc (x, Q 2 )) 2 / (AAf^ 8, (x , Q 2 )) 2 - The optimum set of parameters is found by minimizing x 2 by 
the subroutine MINUIT. 



3 Results 

In our analysis, we use the asymmetry data set of so called "large tables". On the other 
hand, other analyses use small data tables in which Q 2 values are averaged at certain x points. 
Although the present data are not accurate enough to provide Q 2 -dependence information, we 
believe that the raw data should be used as much as we can. For 375 data points with Q 2 > 1 
GeV 2 , we obtain the minimum x 2 values x 2 / d -°- f -= 322 - 6 / 360 and 300.4/360 for LO and NLO, 
respectively. There is a significant x 2 reduction in NLO, so that it is important to analyze the 
data by the NLO expressions. It should be noted that all our NLO results are obtained in the 
MS scheme. There are two major sources for the x 2 reduction, and they are from the HERMES 
proton and E154 neutron data in Figs.|l] and ||. 




Figure 1: Spin asymmetry Ai for the proton. Figure 2: Spin asymmetry Ai for the neutron. 

In these figures, our optimum LO and NLO asymmetry curves are shown at Q 2 =5 GeV 2 . Because 
the data are taken at various Q 2 points depending on the x region, it is not straightforward to 



compare them with the fixed Q 2 curves. However, the figures show that the agreement is 
satisfactory. The errors of deuteron data are so large that its A\ figure is not shown here. 

Obtained polarized parton distributions are shown in Fig||. Because the first moments of 
Au v and Ad v are fixed at positive and negative numbers, respectively, they are mainly positive 
and negative distributions. It is rather difficult to determine the gluon distribution, but a positive 
distribution is favored in both LO and NLO. The antiquark distribution becomes negative, and 
we discuss the issue of its determination in the next section. 

Using the obtained distributions, we show Q 2 dependence of Ai in comparison with data in 
Fig.|j. The distributions are provided at Q 2 =l GeV 2 and they are evolved to larger Q 2 by the 
DGLAP equations. The curves show that there could be strong Q 2 dependence in the small- 
Q 2 region (Q 2 <2 GeV 2 ). Experimentalists sometimes assume that A\ is Q 2 independent for 
evaluating g\ from the A\ data. Although it does not matter due to the present experimental 
accuracy, we should be careful about such an assumption for a precise analysis. 




* Q 2 (GeV 2 ) 

Figure 3: Obtained parton distributions. Figure 4: Q 2 dependence of A\. 

Using these distributions, we get the quark spin content 

AS = 0.201 (LO), 0.051 (NLO). (3) 

The NLO value seems to be significantly smaller than usually-quoted ones 10%~30%. The 
small spin content originates from the small- a; behavior of our antiquark distribution, so that 
we discuss this issue in the next section. 



4 Polarized antiquark distribution 
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Figure 5: Polarized antiquark distributions. 



Because the small- x extrapolation could affect 
the value of the spin content, we compare our 
antiquark distribution with some of the recent 
analyses. In FigJE], our AAC curve is shown to- 
gether with the LSS (Leader-Sidorov-Stamenov) 
and SMC (Spin Muon Collaboration) distribu- 
tions. The antiquark distribution is not explicitly 
calculated in the SMC paper, so that the curve is 
obtained by transforming their distributions. The 
SMC curve deviates from others in the large- x region. However, the large- x difference does not 
matter at this stage because the antiquark distribution does not contribute to g\ at large x. It 
should be clarified for example by the Drell-Yan measurements at RHIC in any case. Another 
important point is in the small- x region. Our distribution falls off rather slowly in comparison 
with others, which results in the small quark spin content. Unfortunately, available data are 
taken in the region x > 0.004, and x is not small enough to determine Aq uniquely. It is also 
still too far away from meaningful flavor decomposition for Aqa 

Due to the lack of small- x data, we had better fix the small- x behavior of Aq. As theoretical 
guidelines, the Regge theory and perturbative QCD could be used. First, the Regge model 



predicts g\{x) ~ x~ a as x — > with the intercepts a of ai(1260), /i(1285), and /i(1420) trajec- 
tories. However, the intercepts are not well known and it is usually assumed as a ai = —0.5 ~ 0. 
The small- x functional form of the GRV98 is xq ~ x~° at Q 2 = 1 GeV 2 according to our 
numerical estimate, so that the Regge prediction is Aq/q — > a; ~ . Second, the perturbative 
QCD could also suggest the small- x behavior. However, we have to assume that singlet and 
gluon distributions are constants as x — > at certain Q 2 (= Q 2 )- Although it is not obvious 
whether such Q\ exists, let us assume Q\ = 0.3 ~ 0.5 GeV 2 . Then, the singlet distribution 
becomes x _012 ~~ ' 09 as x — > 0, hence Aq/q — > x 10 . The Regge and "pQCD" distributions fall 
off much faster than ours at small x. 

From these theoretical suggestions, we ana- 
lyzed the data again by fixing the parameter a.q 
(Aq/q — > x aq ). First, we fixed the parameter at 
1.0 which is suggested by pQCD and it is also 
about the lower bound of the Regge model. Sec- 
ond, it is taken as 1.6 which is the upper bound 
of the Regge. Then, we obtain the NLO results: 




X z =305.8, 
323.5, 



AS 



=0.241, 
0.276, 



for OLn 



=1.0, 
1.6. 



Figure 6: Spin content J AEcte. 



(4) 



The x 2 change for ctq=1.0 from the previous NLO value (300.4) is merely 1.8%, but we notice 
the large change in the spin content (5%— >24%). Because of the small x 2 change, the 0^=1.0 fit 
could be equally taken as a good solution. On the other hand, the x 2 change is rather large for 
a.q=l.Q. According to these results, the spin content is 24%~28% which is in the range of the 
widely-quoted values. In this way, we find that the spin content is very sensitive to the small- x 
behavior of the polarized antiquark distribution and that it cannot be fixed by the present g\ 
data. 



5 Conclusions 

It is clarified in our x 2 analysis of the A\ data that the polarized antiquark distribution is not 
well determined in the small- and large- a; regions. This fact leads to the conclusion that the 
quark spin content is not well constrained only by the present A\ data. We need small- and 
large- x measurements in future. From our studies, we have proposed three sets of distributions: 
LO, NLO-1 with free at q , and NLO-2 with fixed a q = 1.00 
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